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Abstract. This paper explores fast, polynomial time heuristic approxi- 
mate solutions to the NP-hard problem of scheduling jobs on N identical 
machines. The jobs are independent and are allowed to be stopped and 
restarted on another machine at a later time. They have well-dehncd 
deadlines, and relative priorities quantified by non-negative real weights. 
The objective is to find schedules which minimize the total weighted tar- 
diness (TWT) of all jobs. We show how this problem can be mapped into 
quadratic form and present a polynomial time heuristic solution based 
on the Hopfield Neural Network (HNN) approach. It is demonstrated, 
through the results of extensive numerical simulations, that this solution 
outperforms other popular heuristic methods. The proposed heuristic is 
both theoretically and empirically shown to be scalable to large problem 
sizes (over 100 jobs to be scheduled), which makes it applicable to grid 
computing scheduling, arising in fields such as computational biology, 
chemistry and finance. 
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1 Introduction 

With the advent of grid computing, the classical science of scheduling the- 
ory has gained a new sphere of applications. Methods which were originally 
developed for decision making around limited resources in manufacturing and 
service industries have been adopted in the areas of computer science, telecom- 
munication and other computationally intensive disciplines such as computa- 
tional biology, chemistry and finance [5, 22]. Specifically, in the area of com- 
putational finance, where this paper sources its motivation, the problems of 
portfolio selection, pricing and hedging of complex financial instruments re- 
quires an enormous amount of computational resources whose optimal usage 
is of utmost importance to investment banks. The prices and risk sensitiv- 
ity measures of complex portfolios need to be reevaluated daily, for which an 
overnight batch of calculations is scheduled and performed for millions of fi- 
nancial transactions, utilizing thousands of computing nodes. Each job has a 
well-defined priority and required completion time for availability of the re- 
sulting figures to the trading desk, risk managers and regulators. The jobs can 
generally be stopped and resumed at a later point on a different machine which 
is referred to as preemption in scheduling theory. For simplicity of modeling 
the problem, machines are generally assumed to be identical and there is a 
known, constant number of machines available. 

The problem of finding optimal schedules for jobs running on identical ma- 
chines has been extensively studied over the last three decades. Sahni [25] 
presents an O (n log ran) algorithm to construct a feasible schedule, one that 
meets all deadlines, if one exists, for n jobs and m machines. The basic idea 
of the algorithm is to schedule jobs with earliest due dates first, but fill up 
machines with smaller jobs if possible. Note that this method allows the de- 
velopment of an algorithm to compute the minimal amount of unit capacity 
for which a feasible schedule exists. This result has been extended to machines 
with identical functionality but different processing speed, termed uniform 
machines, and jobs with both starting times and deadlines [19]. However, the 
scheduling task becomes more difficult when a feasible schedule does not exist 
and the goal is to minimize some measure of delinquency, often termed tardi- 
ness. Tardiness of an individual job under a given schedule is defined as the 
amount of time by which the job finishes after its prescribed deadline, and is 
considered to be zero if the job finishes on or before the deadline. 

In case of minimizing the maximum tardiness across all jobs, Lawler [14] 
shows that the problem is solvable in polynomial time, even with some prece- 
dence constraints. Martel [19] also used his construction to create a polynomial 
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time algorithm to find the schedule which minimizes maximum lateness. How- 
ever, if our measure concerns the total tardiness instead of the maximal one, 
then even the single machine, total tardiness problem (without weights) was 
proven to be NP-hard by Du et al [7]. A pseudopolynomial algorithm has been 
developed by Lawler [13] for this problem, using dynamic programming, but 
this is for the 1-machine problem and does not have good practical runtime 
characteristics. 

In practical applications, jobs often have relative priorities associated with 
them, represented by positive real weights and the objective becomes min- 
imizing the total weighted tardiness (TWT). Once the NP-hardness of the 
TWT problem was established, most of the research work on the problem 
concerned the development of fast, heuristic algorithms. Dogramaci et al. [6] 
propose a simple heuristic for the total (non-weighted) tardiness problem with- 
out preemption. Rachamadugu et al. [23] then studied the identical machine, 
total weighted tardiness problem without preemption. They proposed a myopic 
heuristic and compared this to earliest due date (EDD), weighted shortest pro- 
cessing time (WSPT) and Montagnes rule on small problem sizes (2 or 5 jobs 
in total) . Azizoglu et al. [3] worked on an algorithm to find optimal schedule for 
the unweighted total tardiness problem without preemption, but their branch 
and bound exponential algorithm is too slow, in practice, for problems with 
more than 15 jobs. Armentano et al. [2] examined the non-weighted problem 
without preemption, and starting from the KPM heuristic of Koulamas [11] 
improved upon it, using tabu search. Guinet [8] applies simulated annealing to 
solve the problem with uniform and identical machines and a lower bound is 
presented in order to evaluate the performance of the proposed method. More 
recently, Sen et al. [26] surveyed the existing heuristic algorithms for the single- 
machine total tardiness and total weighted tardiness problems while Biskup 
et al. [4] did this for the identical machines total tardiness problem and also 
proposed a new heuristic. Akyol et al. [1] provide an excellent recent review 
of artificial neural network based approaches to scheduling problems and pro- 
poses a coupled gradient network to solve the weighted earliness plus tardiness 
problem on multiple machines. The feasibility of the method is illustrated on 
a single 8-job scheduling problem. 

We suggest a novel heuristic for the TWT problem, based on the Hopfield 
Neural Network approach which is shown to perform better than existing 
simple heuristics and has desirable scaling characteristics. Maheswaran et al. 
[16] applied a similar approach to the single machine TWT problem and their 
results were encouraging for a specific 10-job problem. 

In this paper, we first map the problem into quadratic optimization and then 
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the Hopfield net is used to provide fast polynomial time, heuristic solution. 
The topics are organized as follows: in Section 2, we present the problem 
formulation and the model used, in Section 3 existing heuristics are defined 
and explained, in Section 4 our novel HNN approach is introduced, in Section 
5 numerical results are presented and finally in Section 6 some conclusions are 
drawn and directions for future research are outlined. 

2 Problem formulation 

In this section, we give a formal presentation for the problem of optimally 
scheduling jobs on finite number of identical processors under constraints on 
the completion times. The basic formalism is the following: 

• Given N jobs with sizes x = {xi,X2>*3> ■ ■ - x n) £ N N . The processing of 
the jobs can be stopped and resumed at any time, so the processing time 
units of each job need not be contiguous. In the literature this condition 
is known as preemption and also assumes a task started on one machine 
can continue on another [5]. 

• For each job a cutoff time is prescribed by K = {Ki , K2, K3, . . . , K^} £ N N . 
This constraint defines the time within which the job is to be completed. 

• The constant number of processors, the capacity of the system is denoted 
by V G N. 

• We are also given a vector w = {wi, W2, W3, . . . , w^} e M N ,w^ > 0, Vi = 
1 , . . . , N denoting the relative priority (or weight) of each job, which can 
be used in the definition of the objective function. 

A schedule is represented by a binary matrix C € {0, l} NxL where Cy = 1 
if job i is being processed at time slot j, and L denotes the length of the 
schedule. An example is given in (1) where the parameters are the following: 



The first row in (1) denotes the fact that under this schedule C, the first job is 
processed in time steps 1 and 3 (note that preemption is used as the processing 
of this job is not continuous) and therefore the prescribed size 2 of this job 



V = 2, N =3, x = {2,3, 1}, K = {3,3,3}. 
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completes within the prescribed cutoff time of time step 3. Similarly, the 3 
units of the second job complete within the cutoff time of 3 and the third job 
is completed ahead of the cutoff time on time step 2. Summing the columns of 
matrix C, we see that the maximal capacity of V = 2 is fully utilized on each 
time step. 

In order to evaluate the effectiveness of a given schedule C, we define tar- 
diness of a job as follows: 

Ti = max (0, F t - K t ) , (2) 

where F^ is the actual finish time of job i under schedule C: F^ = 
argmax{Ci ; j = 1} (The position of the last 1 in the ith row in scheduling 
i 

matrix C.) 

The problem can now be stated formally as follows: 

N 

C opt := arg min Y~ WiTi. (3) 

c fcr 

Under the following constraints: 

• The sizes of the scheduled jobs in the scheduling matrix are equal to the 
predefined amounts: 

L 

J2Cy=Xi,Vi = 1,...,N. (4) 
J=1 

• The number of scheduled jobs at any given time instant does not exceed 
the capacity of the system: 

N 

^C Lj <V,Vj = l,...,L. (5) 

i=i 

We revisit our previous example with a minor change: V = 2, N =3, 
x = {2, 3, 2}, K = {3, 3, 3} and a weight vector w = {3, 2, 1 }. It can be observed 
that there is no solution, in which all jobs are completed before their cutoff 
times. A minimal weighted tardiness solution is the following: 

/ 1 1 \ 
C = 1 1 1 . 

\ o i o i / 
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3 Existing heuristic methods 

Given the NP-hardness of the scheduling task as shown by Du et al. [7], and 
therefore the amount of time it would take to find the exact, optimal solu- 
tion, in most real- world settings the pragmatic approach of finding a fast, 
sub-optimal, but good solution is followed. As outlined in the introduction, 
this has been a very active field of research, and most studies use the Earliest 
Due Date (EDD) and Weighted Shortest Process Time (WSPT) heuristics as 
benchmarks for evaluating the proposed algorithms. In addition to these, we 
also outline the recently developed Load Balancing Scheduling (LBS) heuris- 
tic and a newly proposed Largest Weighted Process First (LWPF) heuristic. 
Furthermore, we also outline a random processing method which is used as a 
low benchmark for our testing results. 

3.1 Earliest due date (EDD) heuristic 

The EDD heuristic orders the sequence of jobs to be executed from the job with 
the earliest due date to the job with the latest due date. Using the notation 
of Section 2, we relabel the job indices so that the following inequality holds: 

Ki < K 2 < ... < K N . 

Once this ordering is determined, the jobs are allocated to the machines in 
this order, always utilizing the maximum available capacity. Once a job finishes 
and capacity is freed up, the next job using the above ordering is scheduled 
on the freed up machine. It has been shown in [22] that EDD finds the opti- 
mal schedule when one wants to minimize the maximum tardiness on a single 
machine. However, we note that when the objective function includes the rel- 
ative priorities of jobs, EDD is at a severe disadvantage, as it does not include 
consideration of the weights. 

3.2 Weighted shortest processing time (WSPT) heuristic 

The WSPT method is analogous to the EDD method in that it orders the 
jobs according to well-defined criteria and then schedules the jobs according 
to this ordering, utilizing the maximum available capacity available. Once a job 
finishes and capacity is freed up, the next job using the ordering is scheduled. 
Using the notation of Section 2, the jobs are ordered such that the below holds: 

*i_ < X2_ < < x N 
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3.3 Largest weighted process first (LWPF) heuristic 

This heuristic is analogous to the EDD and the WSPT heuristics with the jobs 
ordered according to the following inequality: 

Wi > W2 > . . . > W]\|. 

This is a simple heuristic which allows us get a sense of the importance of 
considering the weights versus the cutoff times. We have not seen this simple 
heuristic mentioned anywhere in the literature, but it turns out to have quite 
good empirical characteristics which is one of the many contributions of this 
paper. 

3.4 Load Balancing Scheduling (LBS) algorithm 

Laszlo et al [12] suggest a novel deterministic, polynomial time algorithm for 
scheduling jobs with cutoff times, in the context of load balancing for wireless 
sensor networks. The basic idea of the algorithm is to start scheduling the jobs 
backwards from the maximal cutoff time, in order of decreasing cutoff times, 
always ensuring full capacity is utilized. For a more formal definition of the 
algorithm, please see [12]. We note that whilst the algorithm has proven to 
be extremely successful for load balancing, it does not take into account the 
weights of the jobs and therefore will suffer from the same shortcomings in our 
setting as EDD. 

3.5 Random method 

In this method, the jobs are ordered randomly and are scheduled in such a way 
as to always fully utilize the maximum available capacity. This unsophisticated 
heuristic is useful as a low benchmark by repeating it a number of times and 
taking the best solution over the multiple repeats. 

4 Novel Hopfield neural network (HNN) based ap- 
proach 

Since the number of possible binary matrices grows exponentially with the 
number of nodes and the length of the schedule, exhaustive search with com- 
plexity (2 N L ) is generally computationally infeasible to solve the optimiza- 
tion problem at hand. Thus, our goal is to develop a polynomial time approxi- 
mate solution by mapping the problem to an analogous quadratic optimization 
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problem, similar to the approach of Levendovszky et al [15] and Treplan et al 
[27]. 

We first review the methods available for optimizing equations in quadratic 
form. Let us assume that matrix W is a symmetric matrix of size n x n and 
vector b is of length n. We seek the optimal n dimensional vector y which 
minimizes the following quadratic function [21]: 

f(y) = -^y T w y + b T y, (6) 

subject to one or more constraints of the form of 

Ay < v, 
By = u. 

If the problem at hand contains only linear constraints then it can be solved 
as presented by Murty et al [20]. In other cases, if the matrix W is positive 
definite, then the function f (y) is convex and the problem can be solved with 
the ellipsoid method presented by Zhi-Quan et al [28]. When W is indefinite, 
the problem is NP-hard (for details see [24]). 

A frequently used powerful, heuristic algorithm to solve quadratic optimiza- 
tion problems is the Hopfield Neural Network (HNN). This neural network is 
described by the following state transition rule: 

/ N , \ 
yt(k + 1) = sgn I Y_ Wij-yj(k) - Si I ,i = mod N k, 

where 

d = -diag (W) , 
W = -W - diag (d) , 
b = b- \d. 

Using the Lyapunov method, Hopfield [10] proved that the HNN converges to 
its fixed-point, as a consequence the HNN minimizes a quadratic Lyapunov 
function: 

1 N N N -. 

c(y) := ~ y. Z ^.jmyj + Y. y& = ~2 yl Wy + fiTy - 

i=1 j=1 i=l 



Thus, the HNN is able to solve combinatorial optimization problems in polyno- 
mial time under special conditions as it has been demonstrated by Mandziuk 



56 



N. Fogarasi, K. Tornai, J. Levendovszky 



[18, 17]. Using this method in practice we have to handle the following prob- 
lem: in some cases the HNN method may get stuck in local minimal point of 
its Lyapunov function. 

These methods motivate us to map the existing optimization problem into 
quadratic form. First the binary scheduling matrix C is mapped into a binary 
column vector y as follows: 

/ Cij C] t 2 ••• Cl,L \ 
02,1 C2,2 • • • C 2 ,L 

\ Cn,1 Cn,2 • • • Cn,i_ / 

y = (Cl,1, Ci^, • • • , Ci ; l, C2,l, • • • , C2,L, Cnj, • • • , Cn^) 1 . 
The original objective function (3) is elaborated as follows: 



mm 
c 



y~ Wi |max ^0, arg max{Cy = 1 } — j 



This objective is equivalent to: 

N 



mm 
c 



1=1 \)=K i+ ] 



The minimization problem is thus equivalent to: 

N L 

C opt := arg min ^_ w ^ C l) ■ 



c 1=1 j=Xt+1 



Therefore the mapping required is: 



2 

The solution is the following: 



1 N L 

^y T W A y + b A T y = Y_ Y_ WtCjj 



W 



A 





1=1 


j=Ki+1 


b A = 


OnLxI 












(: 


D 2 
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where 



D; 



OKjXKj 

°(L-K j )xK j W j ^ ( L— Kj ) x (L— Kj ] 



€ R 



LxL 



Having transformed the objective function to a quadratic form, we now turn 
our attention to doing the same with the two constraints outlined in (4) and 
(5). 

The constraints in (4) can be rewritten as follows: 



3=1 




This will lead to the following equation: 



-y T W B y + b B T y = X Z C ^ 



J=1 



The solution of equation (7) is: 

b B = 2( x llxL x 2lxL 



/ UxL 
UxL 



o \ 


IlxL j 



(7) 



\ o 

Another constraint ensures that the scheduling does not overload the process- 
ing units, described in (5). The required transformation is as follows: 

l 



V 



(8) 



i-l 



i=1 



J=1 



Please note the technicality that this constraint may not be relevant for the 
last few columns where only the remaining jobs have to be scheduled and 
schedule matrices not exhausting the full capacity should not be penalized. 
The total length of scheduling can be written as follows: 

N 

i=l 

V 
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L = max ( L, maxx^, max 

V 1 1 

The number of columns where capacity V needs to be fully utilized is the 
following: 

M = max ^1 , max — f + 1 ^ . (9) 

Taking into consideration (9) the mapping of (8) can be described by the 
following equation: 

1 m / / N 

"2y T W B y + b B T y = Y_ I I Y_ Cy 



V 



t=1 



The solution of this equation is the following: 

be = [Vmx1»0(L-M.)x1j V Mxl,0 (L _ M)xl , . . . , V Mx i,0 (L _ M )xl] , 



/ D D D \ 

D D D 



W 



c 



\ D D D / 

where 

IlVlxM Mx (L-M) 
0(L-M)xM 0(L-M)x(L-M) 
We can combine these mappings into the form of equation (6) as follows: 

,NLxNL 



D 



W = ctW A + (3 W B + yW c € R' 



and 



,NLx1 



b = ab A + (3b B + yb c G R 

Note that the objectives can be controlled with heuristic constants a, (3 and 
y in order to strike an appropriate balance between the weights of different 
requirements. Having this quadratic form at hand, we can apply the HNN to 
provide an approximate solution to the optimization problem in polynomial 
time. 



5 Numerical results 



In this section, the performance of the HNN approach is investigated and is 
compared to the performance of other heuristics. 
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5.1 Simulation method 

Each method outlined in Section 3 and the HNN method outlined in Section 
4 has been tested by simulation on a large and diverse set of input parameters 
with the aim to characterize the algorithms empirically on scheduling problems 
of different size of jobs. The algorithms were implemented in Matlab and tests 
were run in this simulation environment with randomly generated parameters 
such as size of jobs, cutoff times, and weights. The size of each job and its 
cutoff time and corresponding weight are generated as follows: 

xt = random ([1, Ci]) , (10) 

K| = Xi + random ([ci, 1.5 • Ci]) , (11) 

Wi = random ( [1 , C2] ) , (12) 

where random (0) produces a uniformly distributed random integer value in 
range 0. In our simulations the constant Ci equals 10 and C2 equals 5. The 
number of processors (V) is determined as follows: 

V = 0.25 • J. (13) 

The problem is expected to be solved without tardiness when the capacity 
is V = J. Therefore (13) ensures that there is a high likelihood of tardiness 
associated with the generated problem. 

For each problem size, 100 different problems were generated randomly, 
using (10)-(13) and the results the methods were compared in each case. 

The random method was repeated 1000 times for each problem and the 
best solution was used. Furthermore, the HNN method was repeated 1000 
times for each problem with different random starting points and the best 
solution was used. In addition, the heuristical parameters a, (3 and y were 
adjusted between the simulations in order to provide a good balance between 
optimizing the objective function, but also meeting the required constraints 
to produce a valid scheduling matrix. To adjust the heuristical parameters we 
used Algorithm 1. 

Even so, in some cases the HNN is unable to provide a valid solution, because 
the prescribed constraints on job sizes and capacity are violated. Therefore, 
we introduced an error correction function (see Algorithm 2) in order to cut 
and replace the unnecessary 1-s in the scheduling matrix. In our simulations 
parameter e equals 5. All of the results presented in the following sections 
concern valid schedules meeting the required constraints. 
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Algorithm 1 Algorithm for adjusting the heuristical parameters 
Require: x, K, V, e 
<x<-0.1,P ^5,y^5 

repeat 
i <— i + 1 

Ci <-HNN(x,K,V,o,P,y) 

<x^ a + 0.01 
until errors (C^) < e 
for k = 1 -4 t do 

C k <— correct (C k ) 

T k <- calculateTWT (C k ) 
end for 

return min (T) 



Algorithm 2 Correction algorithm for the scheduling matrix produced by the 

HNN 

Require: x, w, K,V, C 
for k = 1 — > L do 

while £{1, C i)k > V do 

Remove 1 from row j, where j is the row in column k with minimal 
weight that has G k = 1 
end while 
end for 

for k = 1 -> N do 

while £[ =1 C k;i > x k do 

Remove 1 from row k from the column I, where I is the righternmost 
column in row k such that C k t = 1 
end while 

while 2I i=1 C k;1 < x k do 

Add 1 to row k in column 1 where I is the lefternmost column where 1 
can be added without violating the capacity constraint V. 
end while 
end for 
return C 
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5.2 Average total weighted tardiness of the different methods 

The first simulation compares the average total weighted tardiness provided 
by the algorithms for different problem sizes. 



t— 




5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 

Number of jobs 



Figure 1: Average TWT produced by each heuristic over randomly generated 
problems depicted as a function of the number of jobs in the problem. 



Figure 1 demonstrates that the best solution achieved by the HNN produces 
better average TWT for all problem sizes than any of the other heuristics. The 
performances of the EDD, LBS and random methods are, on average, worse 
than the WSPT and LWPF heuristics and the HNN for all problem sizes. 
This shows that consideration of the weights in the optimization problem is 
critical to reach near-optimal TWT schedules. It is also clearly visible that the 
HNN consistently outperforms all other methods for all problem sizes. Note 
also that the simple LWPF heuristic, proposed by us, outperforms the WSPT 
heuristic widely used as a benchmark in the literature. 
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5.3 Detailed comparison of the HNN performance versus 
other heuristics 

In this section, we quantify how much better the solution provided by the 
HNN is, in comparison to the EDD, WSPT and LWPF algorithms. Figure 2 
depicts the ratio of average TWT produced by the HNN method versus the 
other heuristics. We conclude that, on average, the best solution of the HNN 
heavily outperforms the traditional solutions: the total weighted tardiness of 
the best HNN solution is 25% - 56% of the EDD, 25% - 84% of the WSPT, 
and 34% — 91% of the LWPF. As the problem size increases, the WSPT and 
LWPF heuristics produce solutions which are closer to the HNN. 



HNN performance vs EDD, WSPT, LWPF 




Number of jobs 

Figure 2: Ratio of average TWT produced by the HNN versus EDD, WSPT 
and LWPF as a function of the problem size 

In order to verify that the HNN consistently outperforms the other heuris- 
tics on a wide spectrum of problems, not just on a few selected problems 
for each given problem size, we ran a more detailed experiment. In this case, 
we investigated 500 randomly generated problems for each problem size and 
computed the percentage of times the HNN produced a better solution than 
LWPF, the next best method. Table 1 shows the result which quite convinc- 
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ingly proves the general applicability of the HNN to different problem types. 
Across the spectrum of problem domain, the ratio is higher than 98.5% for all 
investigated job sizes. 

Table 1: Percentage of problems in which the HNN provides an improved 
solution over the next best heuristic, LWPF 



Job size 


5 10 20 


25 


50 75 


100 




99.9% 100% 99.5% 


99.2% 


99.3% 98.6% 


98.8% 



5.4 Runtime characteristics of the investigated algorithms 

In this section, we summarize the runtime characteristics of the introduced 
algorithms. Table 2 contains the theoretical order of convergence of the algo- 
rithms as a function of the length of the input parameters. 

Table 2: Theoretical order of convergence of the investigated algorithms 



Algorithm 


Order of convergence 


Reference 


Random strategy 


0(L-N] 


[12] 


EDD, WSPT, LWPF 


0(L-N] 


[16] 


LBS 


0(L-N 2 ) 


[12] 


HNN 


L 2 • N 2 ) 


[10], [9] 


Exhaustive search 


O 2 L N ) 





Note that this table shows the theoretical runtime for a single run of the 
specified algorithm, whilst we need to run a constant number of the HNN 
and random iterations. On the other hand, these independent runs can be 
executed in parallel on several computers, multicore machines or even GPU- 
based architectures, so the overall runtime should not be significantly different 
due to the need to run multiple iterations. 

Figure 3 depicts the empirical runtime comparison of the different heuristics 
which confirm the theoretical limits. Although the HNN is the slowest method 
of the investigated heuristics, in practice we see that near optimal solution of a 
1 00 job problem takes only around 6 seconds / iteration (including the iterative 
adjustment of the heuristic parameters as per Algorithm 1) on Core i7@3Ghz 
processor, running non-optimized Matlab code, which is very promising as far 
as its practical applicability in real life applications is concerned. With further 
optimization and parallelization on multicore machines or GPGPU technology, 
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Runtimes 




5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 

Number of jobs 



Figure 3: The average runtime of each algorithm per iteration (including the 
iterative adjustment of the heuristic parameters as per Algorithm 1 for HNN), 
depicted as a function of the number of jobs. Note that the vertical axis is 
using a logarithmic scale 

the HNN method provides better quality solutions than the other heuristics 
within acceptable runtimes. 

6 Conclusions and directions of future research 

In this paper, we studied the NP-hard problem of scheduling jobs with given 
relative priorities (weights) and deadlines on identical machines, minimiz- 
ing the TWT measure. We developed a novel heuristic approach, utilizing 
quadratic programming and the Hopfield neural network and we showed that, 
in general, it outperforms the EDD and WSPT methods. Furthermore, we 
have shown that our approach can be applied in real-time settings for small 
problem sizes and possesses scalability features which are acceptable for real 
world applications. 

More formal methods for the selection of alpha, beta and gamma parameters 
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could be investigated in the future to further improve performance. Another 
idea to explore is whether a more intelligent selection of the initial point for 
the HNN algorithm (eg. the result of the LWPF or WSPT heuristics) would 
improve the algorithms performance over the random starting point which was 
used in our tests. Finally, the performance of the HNN algorithm needs to be 
compared to more complex heuristic methods. 
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